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We report quantum oscillations in underdoped YBa2Cu306.56 over a significantly large range in 
magnetic field extending from w 24 to 101 T, enabling three well-spaced low frequencies at « 440 T, 
532 T, and 620 T to be clearly resolved. We show that a small nodal bilayer coupling that splits 
a nodal pocket into bonding and antibonding orbits yields a sequence of frequencies, Fo — AF, Fo 
and Fo + AF and accompanying beat pattern similar to that observed experimentally, on invoking 
magnetic breakdown tunneling at the nodes. The relative amplitudes of the multiple frequencies 
observed experimentally in quantum oscillation measurements are shown to be reproduced using a 
value of nodal bilayer gap quantitatively consistent with that measured in photoemission experiments 
in the underdoped regime. 

PACS numbers: 71.45.Lr, 71.20.Ps, 71.18.+y 



Discerning the electronic structure of underdoped 
YBa2Cu306+a; in the normal state is a crucial step to un- 
derstanding the origin of unconventional superconductiv- 
ity in these materials [l| . While quantum oscillation mea- 
surements in underdoped YBa2Cu306+a; have revealed 
multiple frequency components [2HJ], ^ has been chal- 
lenging to distinguish an electronic structure from the 
numerous possibilities that can uniquely explain the ob- 
served frequencies. Recent measurements of the chemical 
potential oscillations [f| in underdoped YBa2Cu306+2- in 
strong magnetic fields have helped narrow down these 
possibilities by determining the multiple frequencies to 
arise from a single carrier pocket — likely located at the 
nodal region of the Brillouin zone [|| and supported by 
other experiments sensitive to the density-of-states at 
the Fermi energy 043- An intriguing question there- 
fore arises as to how such a single pocket can give rise to 
the multiple observed frequencies. 

In this paper, we show measurements of quantum 
oscillations in YBa2Cu306+ ;c made using a contact- 
less resistivity technique over an unprecedented range 
in magnetic field extending from rj 24 T to 101 T 
(see Fig. [1]) The large window in inverse mag- 

netic field (A(i) « 0.032 T" 1 ) affords clear resolu- 
tion of three well-separated low frequencies, namely 
« 440 ± 10 T, 532 ± 2 T and 620 ± 10 T — the exper- 
imental limit for distinguishing closely-spaced frequen- 
cies being l/A(-g) « 32 T. We use the unique form 
of the multiple frequency quantum oscillation spectrum, 
in which the central frequency at 532 T is equidistantly 
flanked by two frequencies 532 - 90 T and 532 + 90 T, 
as a clue to infer a possible electronic structure that de- 
scribes the system. We show that an electronic structure 
comprising a single bilayer-split nodal pocket combined 



with magnetic breakdown tunnelling would give rise to 
such a frequency spectrum, and potentially provide an 
explanation for the angular dependence of the observed 
frequencies [H, EH E2 ■ Similar instances of multiple fre- 
quencies arising from magnetic breakdown, are for in- 
stance found in heavy fermion and ferromagnetic families 
of materials 13|, 14 1. A prerequisite for such an expla- 
nation in underdoped YBa2Cu30g + , ! ; is that a reduced 
bilayer splitting at the nodes [l5J persists deep into the 
underdoped regime where the chains are partially occu- 
pied and a form of order 0, Q likely reconstructs the 
Fermi surface. 

The unit cell in YBa2Cu306+a: contains a pair of Cu02 
layers, known as a bilayer (shown Fig. [2ji). Coupling 
within a bilayer [denoted as tj_(k) J sp lits bonding and 
antibonding bands by a finite gap [16j. The general ex- 
pectation, therefore, is for bilayer splitting to transform 
a single frequency Fo of a single layer Fermi surface into 
two separate frequencies describing a bilayer Fermi sur- 
face, where the difference between them is proportional 
to the orbitally-averaged strength of the coupling t± in 
the relevant region of the Brillouin zone. 

The presence of three related low frequencies in the 
experimental data (Fig. [1} comprising two smaller- 
amplitude side frequencies equidistantly spaced from a 
central dominant-amplitude frequency seems difficult to 
reconcile with the simple expectation of two dominant- 
amplitude frequencies from bilayer splitting. An alter- 
native explanation must be sought if the experimentally 
observed sequence of three frequencies is to be explained 
by a single carrier pocket. We find that such a possibility 
arises uniquely at the nodes in the underdoped regime, 
where recent doping-dependent photoemission measure- 
ments show the value of bilayer splitting (;$ 16 meV) to 



2 



N 
X 

CQ 

H 
m 

< 



N 



CO 

c 

CD J> 

■95 o. 



50 

-50 
■100 
■150 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



j — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — r 

_YBa 2 Cu 3 656 « / 




/ \' 










i * 


•* ■> - 

/ V 


5 s * • V : • • : ■ ■ ■ : 

experiment y •• :: • : 

"=2,B MB =14T V V V 


i t J 
I *..' 




rTT 4 > ? =v 







(a) r^Z^f 



30 40 50 60 70 80 90 100 

8 CO 



532 T 




YBa 2 Cu 3 656 b 

experiment 

"=2,S MB =14T • 

--- n=4,B MB =4T . 



440 T 



620 T 



0.25 0.50 0.75 1.00 1.25 1.50 

Frequency (kT) 

FIG. 1: a, Magnetic quantum oscillations (solid line) mea- 
sured in YBa2Cu30e.56 obtained using the contactless re- 
sistivity technique at T — 1.6(1) K consisting primarily of 
w 440 ± 10 T, 532 ± 2 T and 620 ± 10 T frequencies (see b). 
A slowly varying background has been subtracted, the data 
has been multiplied by the monotonic function exp[75T/B] 
for visual clarity. The upper (red) dotted curve is the sim- 
ulated waveform for a nodal pocket split into bonding and 
antibonding orbits of frequency Fo ± AF (using Fq = 530 T 
and AF = 90 T from the experimental data) with n = 2 
magnetic breakdown gaps (see Fig. [3}d) and a characteris- 
tic breakdown field £>mb ~ 14 T. The lower (blue) dotted 
curve is the simulated waveform for a nodal pocket split into 
bonding and antibonding orbits of frequency Fo ± 2AF with 
n — 4 magnetic breakdown gaps (see Fig. 2} and a charac- 
teristic breakdown field Bmb ~ 4 T. Simulations are made 
using Equation ([2]), in which the product R s RoRt is a 
monotonic function of B, and a negative prefactor is used 
for the n = 2 case, b Fourier transform (using a Blackman 
window) of the measured quantum oscillation data compared 
with simulated data. The three experimental low frequencies 
are seen to be reproduced by the simulation with the right 
ratio of amplitudes. 



be greatly reduced compared to that (rj 150 meV) at the 
antinodes [It}. The effect of this significant reduction 
in splitting at the nodes (e g in Fig. [3Jd) in the under- 
doped regime is to introduce the possibility of magnetic 
breakdown tunneling We illustrate this in Fig. [3] 
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FIG. 2: a, A schematic of two YBa2Cu306+:c unit cells 
where the Cu02 bilayers are shaded in pink. Hopping matrix 
elements within (t±) and between (t c ) bilayers are shown, 
b, A contour plot of tx(k) within the Cu02 planes for the 
simplified model expressed as equation |T]), in which case the 
bilayer coupling vanishes at the nodes along which \k x \ — \k y \ 
(shown by dotted gray lines) . 



considering the example of a elliptical nodal hole pocket 
(such as would be created by an ordering vector of form 
Q = (tt, tt)). Bilayer coupling splits the original pocket 
of frequency Fo (Fig. [5^,) into a bonding and antibond- 
ing orbit of frequency Fo — AF and F) + AF respectively 
(Fig. [3Jd) . The effect of the significant reduction in bilayer 
splitting at the nodes, is for a magnetic breakdown orbit 
(dotted line in Fig. to now become observable. As 
seen from Fig.[3Ji, this magnetic breakdown orbit has the 
same frequency (Fo) as the original unsplit Fermi surface. 
Tunneling across a reduced e g leads to the accumulation 
of spectral weight at Fo at the expense of the adjacent 
frequencies Fq — AF and Fo + AF. We thus arrive at 
a possible explanation for experimental observations in 
Fig. [TJ in the case where F) is associated with the cen- 
tral frequency at 532 T and the value of AF is taken to 
be 90 T, thereby yielding additional split frequencies at 
440 T and 620 T. 

The opening of a nodal gap (e g in energy) in 
YE^CusOfl-i-a; found by angle- resolved photoemission 
studies [15j and s upp orted by detailed electronic struc- 
ture calculations [l6|, however, require us to go beyond 
simplified models of the coupling 



<j_(k) = -^-[cos(k x a) — cos(kyb)} 



(1) 



between bilayers which anticipate the splitting be- 
tween bonding and antibonding bands to vanish at the 
nodes 0, [l9[ (occurring along \k x \ = \k y \ in Fig. rjb). 
Possible explanations for the opening of a nodal gap in- 
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FIG. 3: A schematic nodal hole pocket before (a) and after 
(b) bilayer splitting, showing each of the allowed orbits with 
frequencies Fo±AF, and the additional orbit allowed by mag- 
netic breakdown tunnelling (dotted line) with frequency Fo. 
The location of the nodal magnetic breakdown gap e g is shown 
by arrows. 



elude a non vanishing nodal interbilayer couplin g lli 
(t c in Fig. [2^), mixing with the chain bands [It: 
certain forms of Fermi surface reconstruction [2l| 

To ascertain the anticipated relative amplitudes of 
Fo — AF, Fq and Fo + AF frequencies arising from 
nodal bilayer coupling and magnetic breakdown, we make 
use of the Falicov and Stachowiak treatment for mag- 
netic breakdown amplitudes fl8j | . An initial wave of 
unitary amplitude entering a magnetic breakdown junc- 
tion is separated into a transmitted wave with ampli- 
tude A v = i\[P and a reflected wave with amplitude 
A v = y/1 — P. where the tunnelling probability (P) de- 
pends on the magnetic breakdown field -Bmb [25| accord- 
ing to P = e - B " B / B . 

For multiple closed orbits, all possible closed orbits 
are additive, and the contribution from each orbit is 
multiplied by the magnetic breakdown reduction factor 
Fmb = (A/) (A?) v where l v represents the number of 
magnetic breakdown points the orbit traverses by trans- 
mission, and 1^ represents the number of magnetic break- 
down points the orbit traverses by reflection. The form 
of the total quantum oscillation amplitude, therefore, is 
yielded by 



n/2 

N m RMBR s RuRT cos . 

m——n/2 



( 2ir(F Q + mAF) 
B 



(2) 



where N m enumerates closed orbits of the same fre- 
quency Fo + mAF (in the simple case of the elliptical 
nodal pocket considered, n = l v + l n = 2, such that 
m = -1,0,1, and N Q = 2, N ±1 = 1 (table HJ)), while 
i?T, Rd, and R s represent the thermal, Dingle, and spin 
damping factors respectively (26j . 

We find that on introducing magnetic breakdown char- 
acterized by a single parameter Bmb ~ 14 T to a 
simple bilayer-split nodal pocket (shown schematically 
in Fig. EJd), a waveform consisting of three frequencies 



F - AF, F , and F + AF is produced that quan- 
titatively explains two essential aspects of the experi- 
mental data: the field-dependence of the waveform in 
Fig. QJi, and the relative weights of the observed distri- 
bution of frequencies in Fig. Qi>. The simulated data is 
produced using Equation ^ in which the same value 
of -Bmb generates different field-dependent probabilities 
for each of the three frequencies, tabulated in Table |U 
Such a value of magnetic breakdown field corresponds to 
a nodal gap (e g ps \J (eS^MB / tto*) .cf ~ 10 meV (2o|). 
which falls within the upper limit (~ 16 meV) mea- 
sured for the nodal gap in photoemission experiments 
in the underdoped regime [l7|. This is in contrast to the 
significantly larger magnitude of antinodal gap of order 
e g ss 150 meV experimentally measured by photoe- 
mission [l7j . which would preclude magnetic breakdown 
tunnelling, and yield only two split bilayer frequencies as 
discussed earlier. 
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TABLE I: A table counting the number N m of instances 
each orbit m of frequency Fo + mAF and cumulative mag- 
netic breakdown probability -Rmb occurs within the magnetic 
breakdown network. The bilayer-split elliptical (see Fig. [3]b) 
and the diamond (see Fig. [4| pocket scenarios are considered. 

We next turn to the question of whether the reported 
warpins 
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model proposed here. The splitting of the bonding and 
antibonding bands (e a and £b) is given by [lil I20I] 



r(k) a , b = e(k) ± ^/t 2 c +t 2 ± +2t c t 1 _ cos(cfc z 



(3) 



(where c is the c-axis lattice parameter, k z is the in- 
terlayer momentum and e(k) is the single layer dis- 
persion). For comparable values of t c and ij_(k) av- 
eraged over a nodal orbit, equation [3] takes the form 
e(k) a ,b ~ e(k) ± 2tj^coscfc 2 /2 . The modulation along 
k z of this dispersion results in an angular dependence 
of the bilayer-split frequencies p) ± AF which resem- 
bles that of a warped cylinder, potentially consistent with 
the experimentally observed angular dependence of the 
measured frequencies By contrast, the k z modula- 
tion is expected to be strongly suppressed in the case 
of a large disparity in magnitude between t± and t c as 
typically occurs at the antinodal region [lj|. We also 
note that the origin of the three frequencies (Fo — AF, 
Fo and Fo + AF) from magnetic breakdown of a bilayer 
split pocket is consistent with their similar temperature- 
dependence reported in ref. 

Finally, we note that the proposed scheme for mul- 
tiple frequencies arising from a single carrier pocket by 
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FIG. 4: A schematic of a diamond-shaped nodal pocket [a. [24f| 
split by bilayer coupling, yielding 4 magnetic breakdown junc- 
tions with each of the possible orbits shown, (a) shows the 
two orbits (in this case Fo — 2AF and Fo + 2AF) yielded 
by bilayer splitting of the original pocket, (b), (c), and (d) 
show additional allowed orbits due to nodal magnetic break- 
down tunneling, yielding frequencies Fo — 2AF, Fo — AF, Fo, 
Fo + AF and Fo + 2AF (the small breakdown field in the sim- 
ulations of Fig. [T] causes Fo ± 2AF to be weak in amplitude). 



nodal bilayer coupling is generally valid for a pocket at 
the nodal region irrespective of its geometry or specific 
origin. We have thus far considered the simple case of 
elliptical pockets at the node (for instance,_created by a 
spin-density wave or d-density wave [2 , 1||, 22 1 ) with two 
magnetic breakdown points. Another case which can be 
considered (e.g. Fig. 2|) is a nodal pocket with 4 mag- 
netic breakdown points (in the shape of a diamond) , as 
might result from bidirectional charge order 0, [2J] or 
ortho-II band folding [23j]. In this case, the multiplicity 
N m of orbits with frequencies F - 2AF, F - AF, F , 
F +AF, F +2AF is given by 1, 4, 6, 4 and 1 respectively 
(listed in table Q] as seen from orbit counting from Fig.@|, 
corresponding to a modified binomial distribution. For 
this pocket geometry, the waveform and experimentally 
observed ratio of amplitudes of the Fq ± AF and Fo fre- 
quencies is reproduced by a magnetic breakdown field of 
Fmb ~ 4 T [simulation of quantum oscillation amplitude 
after Equation ^ shown in Fig. [1] . The case of a nodal 
pocket of lower symmetry would yield multiple different 
values of AF. 



While previous interpretations of the measured three 
closely-spaced quantum oscillation frequencies attributed 
them to a nodal pocket accompanied by another 
pocket [H, in the scenario presented here for bilayer- 



split frequencies, all three frequencies arise from a nodal 
pocket unaccompanied by another pocket; the latter 
rather than the former scenario is consistent with quan- 
tum oscillation experiments that point to a single carrier 
pocket [f| , in addition to the measured low value of mea- 
sured heat capacity Q. 

In summary, we report quantum oscillation data up 
to 101 T measured on YBa2Cu30g.56, providing a suf- 
ficiently wide range of inverse magnetic field to re- 
veal a clear separation between three low frequencies 
at w 440 T, 532 T, and 620 T. Given the constraint 
that these multiple frequencies must arise from a sin- 
gle carrier pocket 0-0] , here we propose an electronic 
structure involving bilayer splitting at the nodes that 
yields a hierarchy of frequencies similar to experimen- 
tal observations. We show that for a value of mag- 
netic breakdown field (Fmb ~ 10 T) quantitatively con- 
sistent with the photoemission-measured magnitude of 
the nodal gap (e g 16 meV [l7j]) in the underdoped 
regime, a Fermi surface comprising a bilayer-split nodal 
pocket when acted on by magnetic breakdown tunneling 
explains two essential aspects of the data. It explains the 
field-dependent modulation of the waveform and the rel- 
ative amplitudes of the observed sequence of frequencies 
F - AF, F and F + AF. 
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